Arrested Kondo effect and hidden order in URu 2 Si2 
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Complex electronic matter exhibit subtle forms 
of self organization which are almost invisible to 
the available experimental tools, but which have 
dramatic physical consequences. One prominent 
example is provided by the actinide based heavy 
fermion material URu 2 Si 2 . At high temperature, 
the U-5/ electrons in URu2Si2 carry a very large 
entropy. This entropy is released at 17.5 K via 
a second order phase transition jl] to a state 
which remains shrouded in mystery, and which 
was termed a "hidden order" state [2j. Here we 
develop a first principles theoretical method to 
analyze the electronic spectrum of correlated ma- 
terials as a function of the position inside the unit 
cell of the crystal, and use it to identify the low 
energy excitations of the URu2Si2. We identify 
the order parameter of the hidden order state, 
and show that it is intimately connected with 
magnetism. We present first principles results for 
the temperature evolution of the electronic states 
of the material. At temperature below 70 K U-5/ 
electrons undergo a multichannel Kondo effect, 
which is arrested at low temperature by the crys- 
tal field splitting. At lower temperatures, two 
broken symmetry states emerge, characterized by 
a complex order parameter yj. A real yb describes 
the hidden order phase, and an imaginary yb cor- 
responds to the large moment antiferromagnetic 
phase, thus providing a unified picture of the two 
broken symmetry phases, which are realized in 
this material. 

URu 2 Si 2 crystallizes in the body centered tetragonal 
structure displayed in Fig. [l]A Starting form a localized 
point of view for 17-5/ electrons, i.e., treating them as 
core states in density functional theory calculations, one 
obtains a set of wide bands, entirely of spd character, 
as shown in Fig. [Tj3. In an itinerant picture, there are / 
states close to the Fermi level which hybridize with the 
spd bands to form itinerant states of predominantly / 
character, with a bandwidth of the order of 1 eV. This is 
shown in Fig. [Tp. This situation is realized in Ce and Yr 
based heavy fermion materials, whose electronic states 
are well described by narrowing the bands of the density 
functional theory by a factor of 10-1000, to account for 
the heavy mass [3]. 

For URu2Si2, we propose a new scenario for the trans- 
fer of atomic / weight to the itinerant carriers which 
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FIG. 1: (A) Crystal structure of URu2Si2- Gold, grey and 
green spheres correspond to U, Ru, and Si, respectively. The 
black box marks the cut of the unit cell for which we show 
electronic states in Fig. [2] (B,C,D) A sketch of the arrested 
Kondo effect phenomena. In the vertical axis we display the 
energy (in the interval [— 200 meV, 200 me V]) and in the hori- 
zontal axis the momentum (around the center of the Brillouin 
zone). The horizontal dotted line marks the Fermi level. The 
light spd bands are plotted by blue lines, while the bright 
colors correspond to the primarily U-/ spectra. In standard 
band structure calculation, the heavy U bands are either ex- 
cluded from the calculation (shown in B), or included in the 
one electron theory (C). The many body calculation (D) 
shows the new phenomena of Kondo effect developing only 
at intermediate energy scales and intermediate temperatures 
A<T<T*, where A ~ 35K is the crystal field splitting and 
T* ~ 70K is the coherence scale. The temperature in panel 
(D) is T = 19 K< A and color coding of the spectra is the 
projection of the total DMFT electronic spectra to the U-5/ 
character. The DMFT method predicts the U 5/ occupation 
to be 2.0, hence U-5/ configuration in solid is U 4+ . 



we dub the " arrested Kondo effect" . At high temper- 
atures, above the characteristic coherence temperature 
T* ~ 70 K, the U-5/ electrons are localized, and do not 
participate in forming the low energy bands. The U 
atoms settle in 5/ 2 configuration, for which the crystal 
environment chooses the non-degenerate atomic ground 
state. However, the first excited state, which is also non- 
degenerate, is only A = 35 K above the ground state, as 
first observed in polarized neutron scattering experiments 
[3]. Hence in the temperature range A < T < T* , the 
ground state seems doubly degenerate, hence the Kondo 
effect develops, leading to formation of very narrow states 
near the Fermi energy, and a narrow peak in the density 
of states. The Kondo effect is partially arrested below 
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T < A, because the crystal fields splitting between the 
two singlet states induces a partial gaping at the Fermi 
level. The situation is shown in Fig.[lp. A non-dispersing 
slab of / spectral weight is pushed to roughly 8 meV 
away from the Fermi energy While only a small fraction 
of the / spectral weight is present at the Fermi level, 
it has a significant effect, resulting in a mass enhance- 
ment factor which is over 200 at T — 19 K. For temper- 
atures above the crystal field splitting energy, T ~ A, 
the electronic scattering rate is anomalously large. This 
is a signature of multichannel Kondo physics, whereby 
the two singlets play the role of the spin in the Kondo 
problem that scatters degenerate bands. The relevance 
of the two channel Kondo physics to U based impurity 
models was first proposed by D. Cox [5J- Indeed exper- 
imental data on dilute U s Thi_ a; Ru2Si2 systems are well 
described by the two channel Kondo model [BJ. It is re- 
markable that the Fermi liquid regime in URu2Si2 can 
never be reached. Once the two channel Kondo effect is 
arrested by the crystal field splitting, long range order, 
driven by intersite effects, preempts the formation of a 
Fermi liquid state with a low coherent scale. 

The one electron spectral function A(r, to) represents 
the quantum mechanical probability for adding or remov- 
ing an electron with energy to at a point r in real space. 
We calculate it using the combination of the density func- 
tional theory and dynamical mean field theory (DMFT) 

m 

A(t,u) =J2xi*(r) [GwH ~ G^/jaM] Xkfi(r) (1) 

k 

where Xa( r ) are the basis functions, and G a p is the elec- 
tron Green's function expressed in the basis of x- More 
details of the implementation of the DMFT method for 
this problem is given in the online material. 

These computational studies are the theoretical coun- 
terpart of the scanning tunneling spectroscopy technique, 
that has been very useful in describing the properties of 
numerous correlated materials [5] and was very recently 
applied to URu2Si2 [9]. 

The evolution of the electron spectral function A(r,uj) 
along the cut above U atom in the unit cell, as depicted 
in Fig. [TJ\, is displayed in Fig. [2j The curves have an 
asymmetric lineshape of the type 

A {u) oc [(q 2 - 1) + 2q(uj/T)]/[(uj/T) 2 + 1] 

where T is the width, and q measures the asymmetry 
of the lineshape. These curves were first introduced by 
Fano [TU], to describe scattering interference between a 
discrete state and a degenerate continuum of states. The 
continuum of states in scanning tunneling microscopy on 
URu2Si2 is provided by itinerant spd bands, and the nar- 
row discrete states are the U-5/ electronic states, de- 
picted with bright colors in Fig.[Tp. Notice that the line- 
shape has the most characteristic Fano shape when the 
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FIG. 2: (inset) The calculated electronic density A(r, ui = 0) 
in real space at 20 K, as measured by scanning tunneling mi- 
croscopy experiments. We take a representative cut in real 
space, sketched in Fig.l, which cuts the U muffin-tin sphere 
at 3/4 height, guaranteeing a large matrix elements with the 
valence 4/-electrons on [/-atom. The main panel shows the 
energy dependence of the electronic density of states at var- 
ious positions in real space. For each colored dot we show a 
curve in the main panel, which has the same color as the dot. 

position r is on the U or Si atom, with positive asymme- 
try q > on U atom and negative asymmetry q < on Si 
atom. With the exception of the small peak, marked with 
black arrow in Fig. [2] the lineshape on U atom can be 
well fitted to the Fano lineshape with parameter q = 1.24 
and T = 6.82meV. 

The main features of our calculations, including the 
characteristic lineshape, the width, and the strength of 
the asymmetry of the Fano lineshape, as well as the addi- 
tional small peak around 6 meV, were recently measured 
in scanning tunneling experiments by J.C. Seamus Davis 
group [9], the first experiment of this type on a heavy 
fcrmion system. 

To understand the origin of the small peak at 6.8 meV, 
we integrate the spectra over the whole unit cell and re- 
solve it in the different angular momentum channels of 
the different atoms. This orbitally resolved spectra at 
low energies is shown in Fig. [3j\. 

While there are plenty of Si-3p and RuAd state at 
the Fermi level, these states are only weakly energy de- 
pendent in 300 meV interval. On the other hand, the 
heavy U-5/ electron states have a dip and a peak at the 
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FIG. 3: (A) The total density of electronic state (red) and 
the orbitally resolved components at T = 19K. Also shown is 
the U-/ component at higher T = 50 K . Using Fermi liquid 
estimate for the specific hear coefficient 7 oc A(ui = 0)/Z we 
obtain 7 ~ 85 mJ/(mol K 2 ) at T = 50 K (1/Z ~ 60) and 
7 ~ I40mJ/(mol K 2 ) at T = I9K (1/Z ~ 200). (B) The 
DMFT Fermi surface of URu2Si2 in the paramagnetic state 
T = 20K > T or der- We take a two dimensional horizontal cut 
in the the tetragonal body centered Brillouin zone passing 
through the origin. 



same energy of 6.8 meV, which is clearly visible in the 
real space spectra in Fig. [2} Therefore it is the U-5/ non- 
lorentzian shape of the electronic density of state, which 
is responsible for the additional peak and the asymme- 
try variation in the unit cell. This dip-peak lineshape is 
also the key to unravelling the puzzle of the electronic 
structure of the material, since it shows that instead of 
a regular lorentzian Kondo resonance, we have a double 
peak at low temperature, with a pseudogap at zero en- 
ergy. This is the signature of the Kondo effect, loosing its 
strength at temperature T below the crystal field split- 
ting energy T < A . Below the ordering temperature (the 
ordering will be discussed below), the pseudogap in the 
partial U-5/ spectra is considerably enhanced, as pre- 
sented in Fig. [3jV. Note that due to mixing of spd and 
/ states, there is finite admixture of U-5/ state at the 
Fermi level, giving some "heaviness" to the quasiparticle 
even below the ordering temperature. 

The electronic states can also be resolved in momen- 
tum space by computing the spectral function A(k,ui). 
The active state at the Fermi level, given by peaks of 
A(k,uj — 0), determine the Fermi surface of the mate- 
rial, which is displayed in Fig. [3j3. 

For many Ce based heavy fermion compounds, the 
Fermi surface of itinerant DFT calculation, as well as 
DMFT calculation, accounts for the experimentally mea- 
sured de-Hass van Alphen frequencies. On the other 
hand, the DMFT Fermi surface of URu2Si2, shown in 
Fig. [3j3, is qualitatively similar to the localized DFT 
calculation, with 5/'s excluded from the valence bands 
[TT] . Note that the Luttinger's theorem, which counts 
the number of electrons modulo 2, does not constraint 



the Fermi surface of the material. Hence the Fermi sur- 
face does not need to resemble the itinerant DFT result. 

The first principles DMFT calculations demonstrate 
that for URu2Si2 the 5/ 2 configuration has the dominant 
weight. This has important consequences, since it allows 
the physical Fermi surface to have the same volume as 
the system with no / electrons, such as ThRu2Si2 [T2"] . 

Finally, the DMFT Fermi surface is hole-like, and the 
large hole Fermi surface centered at Z displays character- 
istic wave vectors (0.6,0, 0)a* and (1.4, 0, 0)a* shown in 
Fig. [3j3. Recent neutron scattering experiments display 
low energy spectral weight at these incomensurate wave 
vectors [331 OH Q£] . 

In the long range order phase, the Fermi surface sub- 
stantially reconstructs, and multiple small electron and 
hole pockets appear, making this system a compensated 
metal, with the same number of hole and electron carri- 
ers. 

The U-5f 2 configuration has the total angular momen- 
tum J — A and is split into five singlets and two doublets 
in the tetragonal crystal environment of URu2Si2- The 
relative energy of these crystal field levels is still an open 
problem and several sequences have been proposed. For 
a recent review see Ref . [16] . 

The wave functions of the U-5/ 2 configuration, with 
the largest weight in the DMFT density matrix, are 



ID = 



V2 
cos(0) 



(|4>-|-4» 



(|4) + |-4))-sin(^)|0) 



(2) 



(3) 



Here \J Z ) = \J = 4,Jz) is a two particle state of the 
J = 4 multiplct. We obtained </> ~ 0.23-7T. The separation 
between the two singlets is of the order of 35 K. The 
probability for the atomic ground state, |0), and the first 
excited state, |1), at 20K is 0.54 and 0.1, respectively. 

The transition into ordered states requires an under- 
standing of the collective excitations, which are bound 
states of particle hole pairs. The identification of the 
low lying singlets leads us to consider the following order 
parameter for URu2Si2 



i>i = (X 01 (Ri)) 



(4) 



where Xq,i = |0)(1| is the Hubbard operator which mea- 
sures the excitonic mixing between the two lowest lying 
[/-5/ singlets at lattice site Ri. 

This order parameter is complex. Its real part is 
proportional to the hexadecapole operator of A 2g ir- 
reducible representation of the tetragonal symmetry, 
Retjj oc ((J x Jy + J y J x )(J 2 — Jy), introduced in the context 
of nuclear physics long ago. Its imaginary part is propor- 
tional to the magnetization along the z axis, Imip oc (J z ), 
which is the only direction allowed within this crystal 
field set of states. Experimentally, the moment indeed 
points in z direction [?]. 
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At low temperatures, we found two different stable 
DMFT solutions, describing ordered states with non zero 
staggered (ip) with wave vector Q = (0,0,1). The first 
solution has (ip) purely real and describes the hidden or- 
der phase of URu2Si2- This solution has zero magnetic 
moment, does not break time reversal symmetry and has 
nonzero hexadecapole. The second DMFT solution, has 
a purely imaginary (ip) and we associate this phase with 
the large moment antiferromagnet phase, which is ex- 
perimentally realized at pressures larger than 0.7 GPa 
[TTl IT51 Q3E [20]. Hence, the microscopic approach suc- 
ceeds in unifying two very distinct broken symmetry 
states in a single complex order parameter. The existence 
of a solutions with either purely real or purely imaginary 
■0, but without mixtures, indicates a first order phase 
transition between these two phases, as observed in the 
pressure experiments of Ref. [T7J [TSJ HH1 l2"0] . 

It is useful to visualize the meaning of this order pa- 
rameter in a limiting case of a simple atomic wave- 
function, where it takes the form \gs >— cos(6*)|0) + 
s'm(9)e ltp \l) The average magnetic moment of the ground 
state is hence (gs\3\gs) — 4cos(</>) sin(2#) sin(< / 9)*(0,0, 1). 
If tp is 7r/2, we have a phase with large magnetic moment, 
and if cp vanishes, the moment vanishes. 

Our theory provides a natural explanation for a large 
number of experiments, which are very puzzling, when 
examined from other perspectives and suggests new ex- 
periments. 

Firstly, it has been advocated phenomenologically that 
even though the hidden order phase and the large mo- 
ment phase have distinct order parameters, the behavior 
of many observables across the transition is remarkably 
similar. The term adiabatic continuity has been used to 
describe this situation [17j [21], but it is not justified on 
the theoretical grounds since the two phases are sepa- 
rated by a first order phase transition. The proposed or- 
der parameter ipi , which unifies the no-moment and large 
moment phase, explains why even though the two phases 
are separated by a first order phase transition, they are 
in many respects very similar, for example in the critical 
temperature, and entropy change across the transition. 

Secondly, at the hidden order transition, a small gap of 
the order of lOmeV opens in the 5/ quasiparticle spec- 
tra, as seen in optical conductivity [22] . specific heat [23] . 
thermal conductivity |24j . and relaxation rate measure- 
ments [35] [55] . On the other hand the dc resistivity con- 
tinues to decrease at low temperatures [57j, since it is 
dominated by the itinerant spd carriers. This fact is hard 
to understand in a simple itinerant density wave picture. 
The / states gap at low energies is also seen in neutron 
scattering experiments [H [TH [J5] . 

Thirdly, unlike density functional theory calculations 
[28] . our approach described a strongly correlated nor- 
mal state with large entropy, and can account for large 
specific heat coefficient in the paramagnetic state of the 
material. 



It has recently been proposed that one can detect 
hexadecapole order using resonant X-ray technique |30j . 
This would be a direct test of our proposed order pa- 
rameter. A high resolution angle resolved photoemission 
spectroscopy can detect the small kink in the very low 
energy spectra (< lOmeV) of Fig.[lp, which is a unique 
signature of the arrested Kondo effect. It would be in- 
teresting to control the crystal field splitting energy A. 
A slight decrease (increase) of A, will increase (decrease) 
the hidden order transition temperature. 

Finally, our results set the stage for understanding the 
mysterious superconducting transition that takes place at 
the much lower temperatures (T c = 0.8 K). Hidden order 
is fertile ground for superconductivity, while the large 
moment antifcrromagnetic phase completely eliminates 
this instability. Our results suggest that Cooper pairing 
can take place only when the electrons propagate in the 
time reversal symmetric background, but unravelling the 
precise origin of the superconductivity in this material 
will require further sleuthing. 
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ONLINE MATERIAL: ARRESTED KONDO EFFECT AND HIDDEN ORDER IN URU 2 SI 2 

The details of the method 

We carried out the band structure calculation for URu2Si2 using recently developed method which combines the 
Density Functional Theory (DFT) and Dynamical Mean Field Theory (DMFT). The combination DFT+DMFT (for a 
review see Ref. PQ) contains some aspects of band theory, adding a "frequency-dependent local potential" to the Kohn- 
Sham Hamiltonian. It also contains some aspects of quantum chemistry, carrying out an exact local configuration 
interaction procedure by summing all local Feynman diagrams. The details of the current implementation of the 
method are given in Ref. [2] (see in particular Eqs.(2) and (25) that give the site resolved tunneling spectra). 

The crystal structure of URu2Si2 is body-centred tetragonal with point group I4/mmm (ThCr2Si2 type). The 
lattice parameters are a = 4.1240 A, c = 9.5566 A and the internal parameter z — 0.3727 [3]. 

The on-site screened Coulomb repulsion in actinides ranges between 4 — 5 eV. For uranium, which is on itinerant 
side of actindes, our estimate for U is 4.0 eV. The Slater integrals were computed by atomic physics program of R.D. 
Cowan [4], and scaled by 70% to account for screening in the solid. Their values are F 2 = 6.81 eV, F4 = 4.55 eV, 
F 6 = 3.36 eV. 

For the double-counting, we used the standard, localized formula Edc = U(rii mp — 1/2) — J{rii mp /2 — 1/2) with 
J = F 2 /l 1.92196532 = 0.57 eV, and rn mp is the average number of electrons in the U-5/ orbital, for URu 2 Si2 ni mp ~ 2. 

We used two complementary impurity solvers, the One Crossing Approximation (OCA) [TJ |2] and the Continuous 
Time Quantum Monte Carlo (CTQMC) [5J. The first is implemented directly on the real axis and is of invaluable 
help to analytically continue the exact Monte Carlo imaginary time data. For the DFT part, we used the LMTO 
method of S. Savrasov jS] and we crosschecked results with LAPW method, as implemented in Wien2K [7]. 

The URu2Si2 compound is one of the rare cases, where the charge self-consistency of DFT+DMFT is crucial. 
Namely, the DFT potential needs to be computed on the self-consistent DFT+DMFT electronic charge. The non-self 
consistent version of DFT+DMFT severely misplaces the important low energy bands. 
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The broken symmetry states 



At low temperature, we succeeded to stabilize two broken symmetry states with the ordering wave wector Q = 
(0,0, 1): An antiferromagnetic state with large magnetic moment, and a broken symmetry state with no magnetic 
moment. 

There are two independent ways to look at the nature of the broken symmetry state: i) An atomic perspective, 
which was given in the main text of the article. The broken symmetry solution is associated with the local order 
parameter ip = (Xqu) = (|0)(1|), where X^ is the Hubbard operator, which mixes two crystal field singlet states, the 
ground state and the first excited state of the atom. The two atomic states are given by 



10) = ^(|4>-|-4» (5) 
|1) = C ^(|4) + |-4))-sin(0)|O). (6) 



Here atomic states |0) and |4) correspond to the | J = 4, J z = 0) and | J = 4, J z = 4), respectively, ii) Alternatively, 
one can examine the structure of Green's function, the central object of the DFT+DMFT method. To understand 
the symmetry of the broken symmetry solution, one can examine the equal time analog, which is the density matrix 
of the problem. 

Below we will derive the connection between the density matrix of the problem and the atomic eigenstates for the 
case of URu 2 Si2- We first evaluate the projection of the density matrix operator fyfp to the atomic eigenstates. It is 
useful to view the broken symmetry from the point of view of the Green's functions matrix and its equal time limit 
(the density matrix). Exact diagonalization of the atomic problem leads to the following relations 



/5V2/-3/2 


= -o|4)<0| + 6|0)<-4| + -.. 


(7) 


/13/2/5/2 


= -a|0)<4| + 6|-4)<0| + -.. 


(8) 


f 1/2 f -5/2 


= 6|4)(0|-a|0)(-4| + ... 


(9) 


/I5/2/3/2 


= 6|0)<4|-a|-4)<0| + -.. 


(10) 


/I5/2/-5/2 


= Cl |-4)(-4|+c 2 |0)(0| + --- 


(11) 


/I3/2/-3/2 


= C3|-4)<-4|+C4|0)<0| + --- 


(12) 


4/2/3/2 


= C3 |4)(4|+c 4 |0)(0| + --- 


(13) 


f 1/2 f 5/2 


= Cl |4)(4|+c 2 |0)(0| + --- 


(14) 



with a w 0.77, b w 0.25, Ci w 0.97, c 2 w 0.07, c 3 w 0.98, c 4 w 0.63. 

The relation between the crystal field atomic eigenstates (|0), |1), |2)) and the direct atomic states \J,J Z ) (J z — 
4, -4,0) is 



10) 
|1> 
|2> 



I 


— I 


COS <t> 


COS <P 




V2 


sin $ 


sin <1> 


. V2 


V2 





- sin<f> 
cos<I> 



|4) 

1-4) 
10) 



(15) 
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Inverting this matrix equation and inserting it into Eqns. (|7|-( |14[ ), we obtain the followin form of the density matrix 



(/ 5 / 2 /-3/2> 
(/_ 3/2 /5/2> 

(/3/2/-S/2) 
{f -5/2^/2) 

(/-5/2/-5/2> 
(/-3/2/-3/2> 
(4/2/3/2) 
(/s/2/5/2) 



j sin $ 
isin(f> 



(6-0* - <#) H = 



V2 
— i sin $ 

^2~ 
— i sin $ 

— ici cos $ 
2 

— IC3 COS $ 



(a0* — 60) + • • • = 



(a0* 



60). 
(60* - a0) - 
(0*-0) 
(0*-0) 



sin* 

sin* 
sin$ 
sin <f> 



i(6 - a)Re0 + (6 + a)lm0] + • • 
-i(b - a)Rc0 + (6 + a)lm0] + 
i(b — a)Re0 — (6 + a)lm0] + 
— i(6 — a)Re0 — (6 + a)Im-0] - 



103 cos $ 



ici cos <fr 



(0* 
(0* 



0) + 
0) + 



72 

= — ci cos $ Im0 
= — C3 cos $ Im0 
= C3 cos $ Im0 + • • 
= ci cos <I> Im0 + • • 



(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 
(23) 



where is the order parameter = (|0)(1|). Here the dots (■ • • ) stand for the projection to the other atomic states, 
the part of the projection which is nonzero in the paramagnetic state. 

Setting Im0 = and i?e0 ^ 0, we see that S(fl /2 f_ 3/2 ) = -S{f_y 2 fs/2) and <5(/l 5/2 _ 5/2 /_5/2 -5/2) = 
^(4/2 5/2/5/2.5/2)- The symbol S stands for the change of the density matrix compared to its value in the para- 
magnetic state. The structure of the Green's function in the hidden-order phase is therefore 



/ 



A = 





-5/2 


-3/2 


-1/2 


1/2 


3/2 


5/2 


-5/2 


A A 











A £ + A Q 





-3/2 





A B 











A e + A p 


-1/2 








A c 











1/2 











Ac 








3/2 


A £ - A Q 











A B 





5/2 





A e - Ap 











A* 



(24) 



Here A e is the off-diagonal term, present also in the paramagnetic phase and comes due to the tetragonal crystal 
environment. Since the magnetic moment of this state vanishes ((J) =0), we associate this moment-free phase with 
the hidden order state. 

The second solution has Im0 ± and Reip = 0, hence 5{fl /2 f_ 3/2 ) = S{f^ 3/2 f 5/2 ) and S{f_ 5/2 _ 5/2 f_ 5/2 _ 5/2 ) ^ 

^(4/2 5/2/5/2.5/2)- The structure of the Green's function in the local moment antiferromagnetic phase is therefore 
given by 



A 



/ 


-5/2 


-3/2 


-1/2 1/2 3/2 5/2 \ 




-5/2 


A A 





A £ 




-3/2 





A B 


A^ 




-1/2 








Ac 




1/2 








A' c 




3/2 


A e 





A' B 


V 


5/2 





a: 


Q Q A' A ) 



Ull 



(25) 



4 and hence we associate 



with A a 7^ A' A , A B ^ A' B , Ac 7^ A' c . This phase has nonzero magnetic moment 
it with the large moment antiferromagnetic phase in URu2Si2- 

Another way to look at the broken symmetry state is by using the perturbation theory around the atomic limit. 
In the expansion with respect to the hybridization strength A (between the atom and the conduction bands), the 
leading term of the Lutting-Ward functional contains the following terms 



* = X/ C m G mm * [Got * (A_ 5 /2,3/2 - A 3 /2,_5/2) + G10 * (A_ 



•3/2,5/2 - A5/ 2 ,-3/2)l+y^ d m G mm *(G i - G ! io)*(Aj je - A_jJ 
m 

(26) 
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Here the pseudoparticle Green's function Goi and Gio are connected to the Hubbard operators by Goi(r — r') = 
-(T T X i(r)X i(r')) and G i(r - r') = -{T t X 10 {t)X w {t')), respectively. 

One possible self-consistent solution of this functional is Goi = Gio 7^ 0, (A_ 5 / 2i 3/2 — A3/2-5/2) = 
(A_ 3 /2,5/2 — A 5 / 2 ,_3/2) 7^ and Aj z = A_ j3 . This solution corresponds to moment-free phase. 

The second possible solution is G i = -Gio 0, A Jz ^ A- jz and A_ 5/2: 3/2 = A 3 / 2 ,-5/2, A_ 3/2: 5/2 = A 5/2 _ 3 /2, 
which corresponds the antiferromagnetic solution. 



[1] G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Pa rcollet, and C. A. Marianetti Rev. Mod. Phys. 78, 865 (2006). 
[2] Kristjan Haule, Chuck-Hou Yee, Kyoo Kim, arXiv:0907.0195l 

[3] A.P. Ramires, T. Siegrist, T.T.M. Palstra, J.D. Garrett, E. Bruck, A. A. Menovsky, J.A. Mydosh, Phys. Rev. B 44, 5392 
(1991). 

[4] R. D. Cowan, The Theory of Atomic Structure and Spectra (Univ. California Press, Berkeley, 1981). 
[5] Kristjan Haule, Phys. Rev. B 75, 155113 (2007). 
[6] S. Y. Savrasov, Phys. Rev. B 54, 16470 (1996). 

[7] P. Blaha, K. Schwarz, G. K. H. Madsen, K. Kvasnicka, and J. Luitz, Wien2K (Karlheinz Schwarz, Technische Universitat 
Wien, Austria, 2001). 



